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The two-loop box contributions to massive Bhabha scattering may be reduced to two-loop box master integrals 
(Mis) with five, six, and seven internal lines, plus vertices and self energies. The self-energy and vertex Mis may 
be solved analytically by the differential equations (DE) method. This is true for only few of the box masters. 
Here we describe some details of the analytical determination, including constant terms in e — (4 — d)/2, of the 
complicated topology B512m3 (with 5 lines, 2 of them being massive). With the DE approach, three of the four 
coupled masters have been solved in terms of (generalized) standard Harmonic Poly logarithms. 



1. INTRODUCTION 



Bhabha scattering, 



e e 



e+e (7), 



might be used as an extremely accurate luminos- 
ity monitor at the planned International Linear 
Collider ILC, with an accuracy of up to 10 -4 
in the very forward direction. Some more de- 
tails on the general status of this topic may be 
found in |H2l,'ij and in references therein. To meet 
these accuracy needs by theory, one has to control 
the virtual two-loop corrections. This has been 
substantially pushed forward quite recently in 
|4I5| by a derivation of the constant terms of the 
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cross section expansion in the small parameters 
m 2 / s,m 2 /t from a massless calculation [B], com- 
bined with contributions from closed fermion loop 
insertions \7\. Alternatively, we are performing a 
direct and independent calculation of the massive 
two-loop Feynman integrals. Our method relies 
on the calculation of a relatively small number of 
scalar master integrals and the use of algebraic 
relations for expressing all the others. The two- 
loop self energy and vertex master integrals are 
known from [81911011 1] and the references therein; 
for additional informations and an introduction 
to our conventions, see also the webpage [T2] . 
For their determination, the method of differen- 
tial equations has been used |13ll4ll5tTrj) . The 
iterative, nested integration of a system of DEs 
for Mis of rising complexity was possible by use 
of a certain class of functions, the so-called Har- 
monic Polylogarithms (HPLs) ^7] and General- 
ized Harmonic Polylogarithms (GPLs) JSj. This 
works quite nicely for problems with one scale 
(with HPLs) and for some problems with two 
scales (GPLs); see the references mentioned. For 
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two-loop Bhabha scattering, the question arises 
whether a full analytical solution can be obtained 
with HPLs with arguments (+1, 0, —1) and GPLs 
with arguments (+1, 0, — l,y, — 1/y). The answer 
is positive for Mis for self-energies (SE) and ver- 
tices. Presently, we are determining all the box 
masters needed. For the Mis of box type we know 
a few solutions in terms of HPLs and GPLs. How- 
ever, most cases have not been solved so far. One 
reason is the appearance of systems of coupled 
DEs. B613m3 e.g. is a box topology with 6 lines, 
3 of them massive, and one is faced with a sys- 
tem of six coupled DEs. For the topologies B513m 
and B512m3, there are five and four coupled mas- 
ter equations, respectively. 

Here, we present analytical results for topol- 
ogy B512m3, which have been obtained as solu- 
tions of differential equations. B512m3 appears to 
be a complicated and very interesting case. This 
topology is part of decompositions for the topolo- 
gies Bl and B3 |11|. The basic master B714ml 
for the topology Bl is known [T§|, and expressed 
through HPLs and GPLs. So, it is natural to 
assume that its subdiagrams, including B512m3 
should have a similar structure. However, the 
situation might be different. 



2. TWO-LOOP BOX TOPOLOGIES 
AND MASTER INTEGRALS 

There are six double box topologies ■ Three 
of them have master topologies with up to seven 
lines: the planar Bl, the second planar B2, and 
the nonplanar B3. For a reduction, one needs also 
master integrals with six and five internal lines. 
Only few of them are known analytically. The 
completely known double boxes with five lines are 
B514m j^OllO] and B512ml QE]. Further, the di- 
vergent parts of the B512m2 and B512m3 masters 
were determined in ^I] and those of B513m in 
PQ. Here, we will determine the finite parts of 
the three dotted masters for B512m3. 



3. TOPOLOGY B512m3 

The topology B512m3 contributes to the reduc- 
tion of Bl and B3, see Figure It may be ob- 



tained from Bl by shrinking lines 3 and 6, 



B512m3[&i 



d d k l( l d k 2 (fc 2 p 3 )-^ 

D^ 1 D b 2 2 D\ l D\ 5 D b 7 7 



( lJT d/2y 

where the arguments are {b% b 2 64 65 67 bs}- If not 
stated otherwise pi is eliminated by momentum 
conservation, pi + p 2 = P3 + Pi- Further, 

Di = (fci + k 2 + P3 ) 2 - 1, 

D 2 = (fci + k 2 -pi +p 3 ) 2 , 



D 7 



(k 2 
k 2 



Pi -P2+P3) 



1, 



There are four masters, and one may choose those 
shown in Figure [3 

B512m3 = B512m3[111110], 

B512m3dl = B512m3[111210], 

B512m3d2 = B512m3[111120], 

B512m3d3 = B512m3[112110]. 

The singularities of the three dotted masters 
in e = (4 — d)/2 have been determined recently 
using the method of differential equations 
For completeness, we quote them here: 

1 (l + x 2 )y 

B512m3dl = -^Ri^# 

lry(l + x 2 ) 2y 2 (l-x) 2 



el 4(1 



4(1 -yf x{\-y)(l + yf 



(2 



y 



2x(l-y) 2 (l+y) 2 



[2y(l + x 2 ) + x(l - y) 2 ] 

(H[0,y] + 2H[l,y}) 
y(l-x 2 ) 



4x(l - y)' 



;H[0,x] 



y 2 {\-xf 



2x{\-y){\ + yf 



(F[0,0,y] + 2F[0,l,t/])} 



-0(1), 



B512m3d2 



2xy 



xyHlO, x] 



{l-x 2 ){l-yf 



e{l _ x2){l _ y)2 miM-H[-i^x] 

+H[Q 1 Q,x}+ ^+H[0,x}(H[0,y] 
+2H[l,y])) + 0(l), 



3 



3 I 4 

Bl 




Figure 1. The topologies Bl (planar) and B3 
(nonplanar) . 





B512m3 



512m3dl 





B512m3d2 



B512m3d3 



Figure 2. The masters for the topology B512m3. 



B512m3d3 : 



y 



[4C 2 + i/[0,0,y] 
e 2(1 - y z ) 

+2H[O ) l,y]] + 0(l), 



and B512m3 is finite. The functions H are Har- 
monic Polylogarithms (HPLs) in the variables 
x = (V-s + 4 - V^s)/(V-s + 4 + V^s) and y 
(obtained by replacing shy t). A basis for use in 
Mathematica may be found in file HPL4.m at |12) . 

Three of the finite parts of the masters were 
found in several steps. Only B512m3dl and 
B512m3d2 have terms of order 0(l/s 2 ). The 
B512m3d2 decouples in a basis where B512m3dl 
is replaced: {B512m3, B512m3d2, B512m3d3, 
B512m3N}. The B512m3N = B512m3 [11111-1] 



with irreducible numerator is of order 0(1) in 
e. In fact, there are four irreducible numer- 
ators. For the momenta chosen here, they 
are either: (kiP2, fci£>3, ^2^3) and (k±pi or fopi 
or k 2 p2)-, or alternatively: {k\k,2-,k\p%) and 
(fciP3 or k 2 p 3 ) and {k 2 p\ or kip%). Any of 
the numerator integrals may be used to re- 
place a dotted master. The corresponding 
algebraic relation for the chosen case relates 
B512m3N with B512m3, B512m3dl, B512m3d2, 
B512m3d3, Tlllm2 2 , SE312m(s), SE311mONS, 
SE310m(t), V412m2(s), V411m2(t), and will be re- 
produced in 12 . In this basis we may derive 
for B512m3d2[0, x, y] a decoupled first-order s- 
channel DE, which can be easily solved with a 
boundary condition at y = 1: 

—Axy 



B512m3d2[0, X, y] 



{^(2DH[-l/y,x] 



(1 

12H[-y,x] 



x 2 )(l 

- \m\- 



-y) 2 

■l,x] 



-3H[0, x] - UH[0, y] + 2H[1, x] + AH[l, y}) 

+7/4 Cs 

-3H[-l,0,x](H[0,y]+2H[l,y]) 
+2H[0,x](H[l,0,y] + 2H[l,l,y}) 
+(H[0,y]+2H[l,y})(H[-l/y,0,x] 
+H[-y, 0, x] + H[l, 0, x] + H[0, 0, a;]) 
+H[0, 0, y] (H[-l/y, x] - H[-y, x] + H[0, x]) 
+H[1, 1, 0, x] - H[0, 0, 0, y] - 2H[0, 0, 1, y] 
+2H[0, 1, y](H[-l/y, x] - H[-y, x] + 2H [0, a:]) 
+2H[-l/y, -1, 0, a:] - H[-l/y, 0, 0, a;] 
+2H[-y,-l,0,x\- H[-y,0,0,x] 
+H[-1, 0, 0, x] - H[-l, 1, 0, x] - H[0, -1, 0, x] 
+H[0, 1, 0, x] - H[l, -1, 0, x] + H[l, 0, 0, x] 
-3JJ[-l,-l,0,x]}. 

In a next step B512m3d3 can be derived using the 
same basis. However, in the differential equation 
the homogeneous part is absent. This means that 
we are not able to formulate a boundary condition 
for this equation using the analyticity of the so- 
lution. However, an analogous situation appears 
for the differential equation for B512m3d3[0, x, y] 
in the i-channel. So, we have two solutions (which 
must give finally the same numbers), but up to 
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the 'constant' terms c\{y) in the s-channel and 
c>i{x) in the ^-channel. We take x = 1 and get 
in the s-channel a pure function of y. Comparing 
now the two solutions we may determine c\ (y) up 
to an unknown true constant C2(l). For HPLs of 
weight three, this constant is proportional to £3 
and can be fitted numerically: 

B512m3d3[0,x,y] = -±-^{^(4H[-l,y] 

-7/2H [0, y] - 8H[1, y] - 5H[-l/y, x] 
-3H[-y,x] + 4H[0,x]) 
- C 3 /2 

+H[-l,0,0,y] + 2H[-l,0,l,y] 
-5/2#[0, 0, 0, y] - 5H[0, 0, 1, y] - 2H[0, 1, 0, y] 
-4H[0,l,l,y] + tf[-l/y,0,0,z] 
-AH[l,0,l,y]-2H[-l/y,-l,0,x] 
-2H[1, 0, 0, y] + 2H[-y,-l, 0, ar] - H[-y, 0, 0, ar] 
-(H[-l/y, 0, x] -H[-y, 0, x])(H[0, y] + 2H[1, y]) 
-(H[-l/y,x]+ H[-y,x}- H[0,x]) 
{H[0,Q,y]+2H[0,l,y])}. 

Finally, the solution for B512m3dl can be found 
with the original basis {B512m3, B512m3dl, 
B512m3d2, B512m3d3} and using the fact that 
B512m3 is not singular in e (this can be seen by 
IR power counting, but has also been checked 
with a sector decomposition). The vanish- 
ing of the singularity of B512m3 yields an al- 
gebraic relation between B512m3dl[0, x, y] and 
B512m3d3[0, x, y] which may be easily resolved for 
B512m3dl[0, x,y\: 

B512m3dl[0, x, y] = ^i{(3 - 2y + 3y 2 ) 

Zx(i y ) 

H[-l, 0, x] - AyH[0, x] {H[0, y] + 2H[1, y}) 
-(l + yf(H[l,0,x]-H[0,x])} 

+TTn L v iWvO- + + *(1 - y?) 

8x(l - y)-(l +yy 

(H[0, y] - 2H[l, 0, y] + 2H[1, y] - 4JT[1, 1, y}) 
-4(1 - 2x + hx 2 + 2y + Axy + 2x 2 y + y 2 
-2xy 2 + 5x 2 y 2 )H[0,0,x] 
-4(3a; + Ay - Axy + Ax 2 y + xy 2 ) 
(H[0, 0,y] + 2H[0,l,y]) 
-(15 + Ax + llx 2 - Wy + 2Axy 



+ Ux 2 y + Ihy 2 - 2%xy 2 + llx 2 y 2 )( 2 } 

y(l + x 2 ) 
2x(\-y) 2 

TT—^KT v, {2g[-i/y,-i, 0, x] 
-v)v- + vr 

-H[-l/y,0,0,x]+H[-y,0,0,x] 

-2H[-y,-l,Q,x]-H[-l,0,0,v] 
+5H[0,0,l,y]-2H[-l,0,l,y] 
+5/2tf[0,0,0,y] + 2iJ[0,l,0,y] 
+4ff[0,l,l,y] + 2if[l,0,0,y] 
+4ff[l,0,l,y] 
+C3/2 

-(4iJ[0, x] + AH [-1, y] - 3H[-y, x] 
-5H[-l/y, x] - 7/2ff [0, y] - 8H[1, y])C 2 
-(-H[-l/y, x] - H[-y, x] + H[0, x}) 
(#[0,0,y] + 2tf[0,l,y]) 
+ {H[-l/y,0,x]- H [-y,0,x]) 
(H[0,y]+2H[l,y})}. 

The final parts of the dotted masters depend 
on HPLs and additionally on generalized HPLs 
(GPLs) The latter are used here in the no- 
tations of |10I12| . A basis for use in Mathemat- 
lca may be found in file GPL .m at [T3|. While 
B512m3d2[0, x, y] and B512m3d3[0, x, y] have a rel- 
atively simple algebraic structure in x and y and 
depend on functions with definite weights, both 
this is not true for B512m3dl[0, x, y\: this func- 
tion has a complicated algebraic structure and 
depends on HPLs and GPLs of mixed weights. 
However, a comparison with B512m3dl[— 1, x, y], 
the 1/e-term in B512m3dl, shows that both ob- 
servations apply already there. 

Finally, we are left with one unresolved term, 
namely the constant term of the basic master, 
B512m3[0, x, y]. For this, we have built many 
different bases of four Mis, with the aim to 
get finally a solvable differential equation for 
B512m3[0, x, y] with proper denominators (allow- 
ing us to integrate the DE with use of normal 
HPLs and GPLs). However, so far we haven't 
found a basis with a proper structure. It was 
possible to derive two coupled differential equa- 
tions where B512m3[0, x, y] is involved, and then 
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to get a second order differential equation for it. 
We found that this equation has a structure which 
seems not to be integrable with normal HPLs and 
GPLs; the B512m3[0, x, y] might reside in a differ- 
ent class of functions than used so far. 

The masters described here will be posted in 
the Mathematica file Master sBhabha.m at |12| . 
For the calculations we used MATHEMATICA. 

4. SUMMARY 

We have sketched here the derivation of 
analytical solutions for the constant parts 
of Mis B512m3dl[0,x,y], B512m3d2[0,x,y], 
B512m3d3[0,x,y]. The method of DEs was used. 
It is evident that for more complicated cases 
the procedure outlined here will not be applica- 
ble both for principal reasons and for numerical 
ones. E.g., a fitting of constants as it was needed 
for B512m3d3 will in general not be possible due 
to loss of precision of numeric calculations. Fur- 
ther, when the number of coupled Mis is large 
(> 3), the situation starts to be complicated. 
To find solutions, an appropriate basis of Mis 
must be defined and differential equations must 
be prepared by applying differential operators. 
Not always is the choice of Mis with dots and 
without numerators the best one; numerators 
may improve the infrared behavior and change 
the algebraic structure. As mentioned, coupled 
systems of five or six Mis seem to be too compli- 
cated systems to be solved analytically. We don't 
know which class of functions will be finally suffi- 
cient to cover all Mis for Bhabha scattering. The 
topology B512m3 is the most complicated case so 
far, which has been - partly - solved analytically 
with the method of DE. 
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